Upper bounds on r-mode amplitudes from observations of LMXB neutron stars 
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The r-mode oscillations of neutron stars can be potentially powerful probes of cold ultra-dense 
matter, since they couple the low energy degrees of freedom of dense neutron star matter to macro- 
scopic, dynamic observables such as the stellar spin period and temperature. In this paper we present 
upper limits on the amplitude of r-mode oscillations, and their gravitational-radiation-induced spin- 
down rates, in low mass X-ray binary (LMXB) neutron stars under the assumption that the quiescent 
neutron star luminosity is powered by dissipation from a steady-state r-mode. We compute results 
for neutron star models constructed with the APR equation of state for masses of 1.4, 2 and 2.21 
Mq. For the lower mass models (1.4 and 2 Mq) we find dimensionless r-mode amplitudes in the 
range from about 1 x 10 -8 to 1.5 x 10~ 6 . For the accreting millisecond X-ray pulsar (AMXP) sources 
with known quiescent spin-down rates these limits suggest that < 1% of the observed rate can be due 
to an unstable r-mode. Interestingly, the AMXP with the highest amplitude limit, NGC 6640, could 
have an r-mode spin-down rate comparable to the observed, quiescent rate for SAX J1808-3658. 
Thus, quiescent spin-down measurements for this source would be particularly interesting. For all 
the sources considered here our amplitude limits suggest that their gravitational wave signals are 
likely too weak for detection with Advanced LIGO. Our highest mass model (2.21 Mq) can support 
enhanced, direct Urea neutrino emission in the core and thus can have higher r-mode amplitudes. 
Indeed, the inferred r-mode spin-down rates at these higher amplitudes are inconsistent with the 
observed spin-down rates for some of the sources, such as IGR J00291+5934 and XTE J1751-305. 
This can be used to place an upper limit on the masses of these sources if they are made of normal 
nuclear matter, or alternatively it could be used to probe the existence of exotic matter in them if 
their masses were known. 



I. INTRODUCTION 

The r-modes are large-scale currents in neutron stars 
that couple to gravitational radiation and remove energy 
and angular momentum from the star in the form of grav- 
itational waves [IH1] • The physics of these oscillations is 
important in relating the microscopic properties of dense 
matter-such as its viscosity and emissivity-to the macro- 
scopic, observable properties of neutron stars such as 
their spin frequency, temperature, mass and radius. The 
r-modes are unstable to gravitational radiation and their 
amplitudes can grow exponentially if viscous and other 
possible damping mechanisms are not large enough. Be- 
cause they can affect the dynamic properties of the star- 
its spin frequency and temperature evolution- they are 
potentially important probes of the phases of dense mat- 
ter. For example, the boundary of the r-mode instability 
region in the spin frequency - temperature plane, and in 
particular its minimum, which may determine the final 
rotation frequency of the star, is very different for stars 
with different interior compositions [5H7]. 

Whether or not the r-mode instability limits the spin 
rates of some neutron stars it is nevertheless important 
to explore for its potential astrophysically observable sig- 
natures. While mass-radius measurements of neutron 
stars are important for constraining the equation of state 
(EOS) of dense matter, observations of their dynamic 
properties such as spin and thermal evolution are im- 
portant and more efficient in discriminating between dif- 
ferent phases of dense matter. That's because dynamic 
properties are affected by the transport and thermody- 



namic properties of dense matter inside the star, such as 
viscosity, heat conductivity and emissivity which depend 
on low energy degrees of freedom and are very different 
depending on the phase of dense matter present [51 [HI IS]- 
Pulse timing observations of accreting millisecond X- 
ray pulsars (AMXPs) made with NASA's Rossi X- 
ray Timing Explorer (RXTE) are beginning to reveal 
the long-term spin evolution of low mass X-ray binary 
(LMXB) neutron stars. The pulsar recycling hypothesis 
whereby millisecond radio pulsars acquire their fast spins 
via accretion, requires that at least some of these stars 
are spun-up to hundreds of Hz-the current neutron star 
spin record being 716 Hz 1 10] . Still, it remains somewhat 
puzzling that the spin frequency distribution of AMXPs 
appears to cut-off well below the mass-shedding limit of 
essentially all realistic neutron star EOS [TT] , and it has 
been suggested that the r-mode instability may play a 
role in limiting neutron star spin rates |12[ fT5] . Alterna- 
tively, recent work supports the idea that at least some 
of the AMXP population are in or close to "spin equilib- 
rium," essentially determined by the physics of magne- 
tized accretion, and that an additional mechanism, such 
as the r-mode torque, is not required to halt the spin-up 

mug. 

While pulse timing noise in AMXPs has made inter- 
pretation of the X-ray pulsation data difficult, there are 
now several convincing measurements of spin-down dur- 
ing the quiescent phases between accretion outbursts. 
For example, both SAX J1808.4-3658 (hereafter SAX 
J1808) and IGR J00291+5934 (hereafter IGR J00291) 
show spin-down of the neutron star between outbursts, 
at rates, i>, of approximately 1 — 3x 10~ 15 Hz s _1 |161I17| . 
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The most convincing evidence for an accretion-induced 
spin-up during outbursts is in IGR J00291, for which a 
peak value of about 3 x 1CP 13 Hz s~ x has now been in- 
ferred [16 . The magnitude of the "instantaneous" spin-up 
torque is larger in this source than the spin-down torque 
by about a factor of 100. The long-term evolution is then 
a competition between the accretion outbursts which spin 
it up and the quiescent spin-down intervals. If the out- 
bursts are frequent enough then the star will be spun-up, 
if not, then long-term spin-down occurs. As the spin- 
up torque is proportional to the mass accretion rate, the 
outcome can also be expressed in terms of the long-term, 
average mass accretion rate. 

The observed quiescent spin-down rate, v sc i, puts an 
upper limit on the torque that can be present due to any 
r-mode, as it has to be less than the observed torque 
of 27r/!> s( 2, where / is the moment of inertia of the neu- 
tron star. Now, spin-down of pulsars is typically ascribed 
to the magnetic-dipole torque that is almost certainly 
present at some level due to the large scale magnetic field 
of the neutron star. Indeed, magnetic field estimates are 
usually obtained by equating the entire observed spin- 
down to that expected theoretically for magnetic-dipole 
radiation. But the observed spin-down is due to the sum 
total of any torques present. Assuming that both r-mode 
and magnetic-dipole torques are present, then the magni- 
tude of the r-mode torque is equal to the observed torque 
minus the magnetic-dipole torque. Since the magnetic 
field strength and geometry of AMXPs are typically not 
known independently from spin-down estimates, the pre- 
cise value of the magnetic-dipole torque is not known a 
priori. However, to the extent that the magnetic torque 
accounts for the majority of the observed spin-down, as 
is typically assumed, then the r-mode torque is likely to 
be much less than the measured spin-down torque. 

In addition to torquing the neutron star the viscous 
damping of the r-modes acts as an internal source of heat. 
Calculations of the coupled thermal and spin evolution 
of accreting neutron stars including the effects of r-mode 
heating and gravitational radiation have been carried out 
by several authors [7J [TBI [T5] . The primary assump- 
tion in these analyses has been that the long-term accre- 
tion (spin-up) torque is balanced by the r-mode torque 
due to the emission of gravitational radiation that carries 
away angular momentum. That is, the long-term aver- 
age accretion torque is in equilibrium with the r-mode 
torque. As noted above, the recent pulse timing observa- 
tions of AMXPs indicate that this equilibrium assump- 
tion is likely not realized in practice. Moreover, Brown & 
Ushomirsky (TSJ showed that if the average, fiducial ac- 
cretion torque given by N acc =< M > (GA/i?) 1 / 2 , where 
< M > is the long-term average mass accretion rate, were 
balanced by the r-mode torque then the quiescent lumi- 
nosities of some accreting neutron star transients should 
be substantially larger than observed due to heating from 
the r-modes. Here we use a similar argument to place up- 
per bounds on the r-mode amplitudes that can be present 
in accreting neutron stars assuming that r-mode heating 



provides the source of neutron star luminosity in the ab- 
sence of accretion. We then use these amplitude limits to 
assess the level of r-mode spin-down that can be present, 
its relation to observed spin-down rates when available, 
and the expected strength of gravitational radiation. 

The paper is organized as follows. In §2 we review the 
basic theory of the r-mode instability and how it couples 
to the thermal and spin evolution of neutron stars. In 
§3 we describe our methods for constraining the r-mode 
amplitudes using the observed properties of quiescent, 
LMXB neutron stars, and we assess the implications for 
r-mode spin-down and the emission of gravitational ra- 
diation. We provide a summary and discussion of the 
implications of our findings in §4. 

II. SPIN-DOWN DUE TO UNSTABLE R-MODES 

The "flow pattern" of the r-modes is prograde in the in- 
ertial frame and retrograde in the rotating frame, which 
means that it moves in the same direction as the star's 
rotation as seen by an observer at infinity, but in the op- 
posite direction as seen by an observer at rest on the star. 
Any mode that is retrograde in the co-rotating frame 
and prograde in the inertial frame grows as a result of 
its emitting gravitational waves. This is the well-known 
Chandrasekhar-Friedman-Schutz (CFS) mechanism, and 
it means that gravitational radiation drives the mode 
rather than damps it. Viscosity on the other hand tends 
to damp the r-mode by transferring angular momentum 
from the mode to the rigid (unperturbed) star. The total 
angular momentum of a perturbed star can be written as: 

j = m-j c (i) 

where I is the moment of inertia and Q is the angular 
rotation frequency of the star and J c is the canonical 
angular momentum of the mode, given by: 

J c = -^JMR 2 Qa 2 (2) 

where a is the amplitude of the r-mode and J is a 
dimensionless constant defined by |20| : 

where p is the run of density within the neutron star. 
The moment of inertia of the star, /, can also be written 
as / = IM R 2 where / is defined by: 

is w^C"' Ad '-- (4) 

To derive the equations for the dynamical evolution of 
the star we use the following argument by Ho & Lai [2T] : 

The canonical angular momentum of the mode in- 
creases through gravitational radiation and decreases by 
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transferring angular momentum to the star through vis- 
cosity 



dJc 
dt 



-Jc 



TV 



■Jc 



(5) 



where the viscous damping time Ty is given by — = 

— + — + .... Here, t<? and tr refer to shear and bulk vis- 
cosities, respectively, and the ellipsis denotes other pos- 
sible dissipative mechanisms, such as boundary layer ef- 
fects |351|55]. 

The second evolution equation is obtained by writing 
the conservation of the total angular momentum J of the 
perturbed star which says that the total angular momen- 
tum of the star decreases due to gravitational radiation 
and increases due to accretion: 



dJ 

lit 



2Jc 

T~G 



(6) 



where N acc is the accretion torque. For a fiducial 
torque iV acc can be written as M(GMR) 1 / 2 which as- 
sumes that each accreted particle transfers to the star 
an angular momentum equal to the Keplerian value at 
the stellar radius R [TB]. In the previous equations the 
quantities — are given by, 



masses of 1.4M© and 2.OM only allow modified Urea 
neutrino emission in the core, but the one with a mass of 
2.21M allows direct Urea neutrino emission in a core of 
radius 5.9 km. Direct Urea processes are very sensitive 
to the proton fraction of dense matter. The required pro- 
ton fraction is roughly 14% in the case of the APR equa- 
tion of state, reached at relatively high density n ~ 5n , 
where no is the nuclear saturation density, this could be 
different for other equations of state [S]. 

The evolution equations for the amplitude of the r- 
mode, a, and spin frequency of the star, fi, can be writ- 
ten by using equations [5] and [6] and substituting J c from 
Eq. [2j The third equation -which describes the tempera- 
ture evolution- can be obtained by noting that the tem- 
perature of the star decreases due to thermal emission 
from the surface and neutrino emission from the interior, 
which in an average mass hadronic star is dominated by 
modified Urea processes (in a massive star direct Urea 
processes can also happen in the core and should be in- 
cluded in the neutrino emissivity as well) , and it increases 
due to the viscous dissipation of the r-mode energy, Py. 
This gives the following equations for the evolution of 
spin frequency, r-mode amplitude and temperature, 



1 



2E r 



(7) 



where E c is the canonical energy of the r-mode, Pq is 
the power radiated by gravitational waves and Pb and P$ 
are the powers dissipated due to bulk and shear viscosity, 
respectively (in natural units where c = h = fee = 1)) 



Pg 

Pb 
Ps 



1 



a 2 n 2 JMR 2 



(8) 
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16to 
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(2m + 3)(m + 1) 5 k 2 
(m - l)(2m + l)S m A 3 (+ c a D R 3 a 2 n 2 



(9) 
(10) 

(11) 



Here we consider m = 2 r-modes, and Aqcd and Aew 
are characteristic strong and electroweak scales intro- 
duced to make V and S dimensionless. In our calcula- 
tions we have used Aqcd = 1 GeV and Aew — 100 GeV. 
The dimensionless parameters V , S, I and J, which in- 
volve radial integration over the star, and 5 and a are 
given in Table |TJ for three different neutron star models 
that we study in this paper. All three models here are 
made of normal hadronic matter with the APR equa- 
tion of state [21] but they have different masses (1.4 M , 
2.0 M Q and 2.21 M ) and radii [5]. The two models with 
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dE 
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(12a) 
(12b) 
(12c) 



where Q = |j, and the viscous dissipated power is 
Py = Ps + Pb + . . • , and again the ellipsis denotes other 
possible dissipative processes. 

Here, L v . L 7 and Cy, which are the total neutrino 
luminosity, photon luminosity and specific heat of the 
star, are given by (in natural units) 



JlGM 2 R 2m+2 a 2 n 2m + 4 



L„ 



AttR^A 



3 A 9—9 
QCD 



L 



A 4 



rpV 



L 1 = 47ri? 2 ( 7T 4 // 



C v = 4ttA q - c v d R 3 C v T v , 



(13) 

(14) 
(15) 



where T is the core temperature of the star and T e f f 
is the surface temperature. The dimensionless parame- 
ters L v and Cy (that involve radial integration over the 
star) and and v are given in Table [I] for different stellar 
models. 

Fig. [T] shows the r-mode instability window computed 
for a 1.4Af neutron star using the APR EOS (the modi- 
fication to the instability window is rather modest in the 
case of the two other stellar models considered in this pa- 
per [5]). In the region above the curve the gravitational 
radiation timescale is smaller than the viscous damping 
timescales, therefore r-modes are unstable and their am- 
plitudes grow exponentially in this region. As can be 
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Neutron Star 
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n K [Hz] 
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NS 1.4 M Q 


core 
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1.81 x 10~ 2 
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2.36 x 10 -2 
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6 
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0.300 


2.05 x 10" 2 
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1.16 x 10" 3 
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NS 2.21 M Q 
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2.02 x 10 -2 


5.05 x 10" 4 


9.34 x 10 -4 


2.62 x 10 -2 


1.29 x 10 -2 
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Table I: Radius, Kepler frequency and radial integral parameters that appear in the moment of inertia, angular momentum 
of the mode, dissipative powers due to shear viscosity (from leptonic interactions) and bulk viscosity (due to Urea processes) , 
specific heat and neutrino luminosity for different neutron star models considered in this work [5]- 



seen in this figure most of the LMXB's are in the unsta- 
ble region for normal hadronic stars where the damping 
is due to the shear and bulk viscosity from leptonic scat- 
tering and Urea processes, therefore there must be some 
non- linear mechanism that saturates the amplitude of the 
unstable r-modes at a finite value. Supra-thermal bulk 
viscosity (25] is one of these non-linear mechanisms, but 
in the case that only the core of the star is considered and 
the effects of the crust are ignored, it can only saturate 
the r-mode at large amplitudes |26| . Mode coupling is an- 
other mechanism that can saturate the r-mode at small 
amplitude when hyperons and therefore large hyperon 
bulk viscosity is included in the computations |27l |2"5] , 
At this point it is not entirely clear which mechanism is 
actually responsible for saturating the r-mode amplitude, 
however, in this paper our primary interest is to obtain 
upper bounds on r-mode amplitudes from observations 
of neutron star transients and this does not require a 
precise understanding of the saturation mechanism. Un- 
derstanding the detailed physics of the saturation mech- 
anism is an important issue, but it is beyond the scope 
of this paper. 

Here we note that if we were to consider the existence 
of exotic forms of matter in the star, then the shape of 
the instability window would change due to their differ- 
ent shear and bulk viscosities |5J [7J The shape of 
the instability window may also be different due to crust 
boundary layer effects |Bj, but even in that case for realis- 
tic models most of the LMXB sources are in the unstable 
region. Therefore, we assume that these sources have 
unstable r-modes which are emitting gravitational waves 
but the amplitude of the r-mode is not growing exponen- 
tially, and there is a mechanism that can saturate the 
growth of the r-mode. In addition, we assume that all 
of these stars are made of normal hadronic matter (con- 
structed with the APR EOS) and all of the sources that 
are in the unstable region for hadronic stars in Fig. [T] are 
emitting gravitational waves at a constant amplitude. 

Since we don't know which mechanism is actually re- 
sponsible for saturating the r-mode amplitude, we simply 
assume that there is a nonlinear mechanism that satu- 
rates the mode. When a hits the saturation amplitude, 
the right hand side of Eq. |12b| becomes zero which im- 
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Figure 1: The r-mode instability region for a 1.4Mq NS 
constructed with the APR equation of state in the spin fre- 
quency vs. core temperature plane. Also shown are some of 
the low mass X-ray binaries which have been considered in 
this work. The horizontal line extending rightward from the 
temperature symbols (the black squares) shows the difference 
between two models for relating the surface temperature to 
the core temperature (i.e. the difference between a fully or 
partially accreted envelope). The difference between the core 
temperatures in these two cases gets larger as the surface tem- 
perature increases, but even for the highest surface tempera- 
ture considered in this work the difference is not big enough 
to change our results. 
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since 



where we have neglected the last term in Eq. 
it is much smaller than the other terms. 

Therefore, when the amplitude is saturated, in all 
of the evolution equations ^ should be replaced by 

^( x=^2q ) and \ p v\ by ^.g . Since viscosity alone 
can not stop the growth of the r-mode and an extra 
mechanism is needed to do that, therefore at saturation 
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— (i — T7=r) is larger than — and the reheating term on 
the right hand side of the temperature evolution equa- 
tion (Eq 12c I will be larger than the reheating due to 



bulk and shear viscosity jS] 



Source 


A = t 


(Hz s- 1 ) 


Osp.eq 

(1.4M ) 


(2.OM ) 


a.p.e, 
(2.21M ) 


IGR J00291 


13 
1363 


5 x 10~ 13 


1.46 x 10~ 7 


1.22 x 1CT 7 


i.4i x itr 7 


SAX J1808-3658 


40 
2x365 


2.5 x ltr 14 


3.20 x 1(T 7 


2.66 x l(r 7 


3.08 x ur 7 


XTE J1814-338 


50 
19 X 365 


1.5 x 10" 14 


2.12 x itr 7 


1.76 x 1(T 7 


2.04 x 10~ 7 



Table II: Duty cycle, spin-up rate |30| and upper limits on 
r-mode amplitudes from the "spin-equilibrium" condition for 
three neutron star transients with assumed masses of IAMq, 
2.0M© and 2.21M©. 



III. CONSTRAINING THE R-MODE 
AMPLITUDE 

A. Constraints from "spin equilibrium" 

Here we compare two different methods for constrain- 
ing the r-mode amplitude, a, from the observations of 
LMXB neutron star transients. The first one, which 
gives larger values for a, is based on the spin equilib- 
rium assumption O (T5J HH1 E2] where we assume that in 
an outburst-quiescence cycle all the spin-up torque due to 
accretion during the outburst is balanced by the r-mode 
spin-down torque due to gravitational radiation in the 
whole cycle. This is similar to the prescription consid- 
ered by previous authors (0HB])) but rather than using 
a "fiducial" torque estimated from the long-term average 
M, we can now use the observed spin-up rates and out- 
burst properties to directly constrain the torque. 

Therefore we have: 



2-kIvA = 



2Je 



(17) 



where v is the spin-up rate during outburst and A = t ^ L 
is the ratio of the outburst duration, t 0l to the recurrence 
time, t r . We can estimate the left hand side of Eq. [17] 
from X-ray observations of LMXBs and since the right 
hand side is a function of a through J c (see Eq. |2| we can 
determine a. In Table [TT| we give the results for a for three 
sources for which there are now estimates of the spin- 
up rate due to accretion. The values of a are given for 
the three different neutron star models considered in this 
work, and they are all in the range from 1 — 3 x 10~ 7 . 

At these amplitudes the inferred r-mode spin-down 
rate would be competitive with the magnetic dipole spin- 
down rate which almost certainly exists in these LMXBs, 
and which is quite likely the dominant spin-down mech- 
anism. Moreover, the amplitudes are also comparable 
to those deduced assuming r-mode spin equilibrium with 
the "fiducial" accretion torques estimated by Brown & 
Ushomirsky [TS]. These authors also demonstrated that 



at such amplitudes some of these objects should have 
significantly higher quiescent luminosities due to r-mode 
reheating than observed. These results suggest that the 
r-mode torque does not balance the accretion-driven spin- 
up torque and that r-mode amplitude estimates based on 
the "spin equilibrium" assumption will overestimate the 
true amplitude. 



B. Constraints from "thermal equilibrium" 

Here, we use the same thermal equilibrium argument 
outlined by Brown & Ushomirsky, but rather than esti- 
mating the quiescent luminosity using the r-mode ampli- 
tude deduced from "spin equilibrium," we use observa- 
tions of the quiescent luminosity of LMXBs to directly 
constrain the amplitude of the r-mode. This works be- 
cause in a steady-state, gravitational radiation pumps en- 
ergy into the r-mode at a rate given by, Wd — (1/3)0 J c = 
—2E c /tg- This expression has the familiar relationship 
for a power dissipated by an applied torque, and in this 
case it is simply the r-mode torque due to gravitational 
radiation. In a thermal steady state all of this energy 
must be dissipated in the star. Some fraction of this heat 
will be lost from the star due to neutrino emission and 
the rest will be radiated at the surface. We further as- 
sume that all of the heat emitted from the surface of the 
star during quiescence is due to the r-mode dissipation in- 
side the star. The resulting r-mode amplitude limits are 
upper bounds in the sense that the observed luminosity 
reflects the contribution from r-mode heating as well as 
any additional sources of heat that are present, such as 
for example due to accretion and the nuclear processing of 
accreted material. If any such sources of heat are present 
then the actual r-mode amplitude will be less than the 
upper bounds given here. For the sources that we study 
in this work, since we know the values of surface tem- 
peratures and quiescent luminosities from observations, 
we can estimate the core temperature and therefore de- 
termine the neutrino luminosities to estimate the total 
amount of heat deposited in the core of these systems by 
gravitational radiation. 

To compute the core temperatures we use equation A8 
in [33] which relates the effective surface temperature of 
the star T" e // to the internal temperature T&, which is the 
temperature at a fiducial boundary at pf, = 10 10 g/cm 3 
for a fully accreted envelope and is valid at not too high 
temperatures (Tf, < lO 8 ^), 



( T eff u 



10 



14, 



r)(18.1 



10 9 K 



(18) 



where g = GM/(R 2 ^Jl — r g /R) is the surface gravity 
and r g = 2GM/c 2 . Here we assume that the neutron 
star's core is isothermal and since the thermal conduc- 
tivity of the crust is high |36| we have T core = Th to 
good approximation. If we instead use equation A9 in 
[35 which gives Tf, for a partially accreted envelope, and 
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Source 


v a (Hz) 


Distance 


kT P t f 






L 1 {erg.s 1 ) 


L v (erg.s 1 ) 


L 1 (erg.s 1 ) 


L v (erg.s ) 






(kpc) 


(eV) 


(1.4M a ) 


(2.21M Q ) 


(1.4M©) 


(1.4M Q ) 


(2.21M ) 


(2.21M Q ) 


4U 1608-522 


620 


4.1 ±0.4 


170 


1.25 x 10 s 


1.34 x 10 s 


1.21 x 10 34 


2.16 x 10 32 


3.03 x 10 34 


4.52 x 10 38 


IGR J00291+5934 


599 


5± 1 


71 


2.96 x 10 7 


3.18 x 10 7 


3.65 x 10 32 


2.06 x 10 27 


9.27 x 10 32 


7.94 x 10 34 


MXB 1659-29 


556 


~ 10 


55 


1.96 x 10 7 


2.07 x 10 7 


1.35 x 10 32 


7.64 x 10 25 


3.29 x 10 32 


6.08 x 10 33 


Aql X-l 


550 


4.55 ± 1.35 


94 


4.70 x 10 7 


5.06 x 10 7 


1.12 x 10 33 


8.40 x 10 28 


2.87 x 10 33 


1.31 x 10 36 


KS 1731-260 


524 


7.2 ± 1.0 


70 


2.89 x 10 7 


3.12 x 10 7 


3.44 x 10 32 


1.70 x 10 27 


8.88 x 10 32 


7.13 x 10 34 


XTE J1751-305 


435 


9±3 


< 71 


2.96 x 10 7 


3.18 x 10 7 


3.65 x 10 32 


2.06 x 10 27 


9.27 x 10 32 


7.94 x 10 34 


SAX J 1808-3658 


401 


3.5 ±0.1 


< 30 


7.23 x 10 6 


7.69 x 10 6 


1.21 x 10 31 


2.63 x 10 22 


2.99 x 10 31 


1.60 x 10 31 


XTE J1814-338 


314 


6.7 ±2.9 


< 69 


2.82 x 10 7 


3.01 x 10 7 


3.24 x 10 32 


1.39 x 10 27 


8.12 x 10 32 


5.72 x 10 34 


NGC 6440 


205 


8.5 ±0.4 


87 


4.11 x 10 7 


4.46 x 10 7 


8.10 x 10 32 


2.88 x 10 28 


2.11 x 10 33 


6.10 x 10 35 


XTE J1807-294 


191 


8.35 ±3.65 


< 51 


1.72 x 10 7 


1.83 x 10 7 


9.84 x 10 31 


2.71 x 10 25 


2.46 x 10 32 


2.95 x 10 33 


XTE J0929-314 


185 


7.8 ±4.2 


< 50 


1.66 x 10 7 


1.79 x 10 7 


9.06 x 10 31 


2.06 x 10 25 


2.31 x 10 32 


2.53 x 10 33 



Table III: Spin frequency, distance to the source |31| . effective temperature at the surface of the star |32H34| . core temperature, 
photon luminosity at the surface of the star and neutrino luminosity for both IAMq & 2.2M© NSs. Note that kT e ff given 
in this table is for a IAMq NS with a radius of 10 km but in computing the core temperatures and luminosities for different 
neutron star models the appropriate red-shifts have been used. We note that for the sake of brevity the core temperatures for 
the 2Mq models are not included in the table, however, their values are all « 5% less than those for the 2.21M© model. 



a column depth of ^ = 10 9 g/cm 2 [7] we get core tem- 
peratures slightly higher than those derived from Eq. [18] 
The right hand side of the error bars on the temperatures 
in Fig [I] shows this difference. It is really only relevant 
for a single source, 4U 1608-522, but even in this case it 
is less than a 50% increase, and the difference is always 
small enough that it doesn't qualitatively change our re- 
sults in the remainder of the paper. To compute Tf, we 
have used the effective surface temperatures, T e ff, given 
by [32H53] 1 . Note that those temperatures are computed 
for a IAMq NS with a 10km radius and we have used 
the appropriate redshifts to compute T core for our neu- 
tron star models. 

Having surface and core temperatures for these sources 
we can use Eqs. [T3]and[l~4]to evaluate their neutrino and 
thermal luminosities for different stellar models. The val- 
ues of the core temperatures, photon and neutrino lumi- 
nosities for different sources for 1.4M and 2.21M NS 
models are computed and given in Table |III| Note that 
in the case of 1.4M NSs the neutrino luminosity is only 
due to the modified Urea reactions, but for the 2.21M 
NS it is due to both direct and modified Urea reactions. 
By comparing the neutrino and thermal luminosities in 
Table |III| one can see that at temperatures relevant for 
the LMXBs the standard neutrino cooling from modi- 



1 It should be emphasized that the temperatures given in Table 2 of 
Heinke et. al. 2007 & 2009 32, 33] are effective temperatures at 
the surface of the star and not redshifted surface temperatures 
seen by an observer at infinity. We note that they have been 
incorrectly assumed to be redshifted temperatures in references 

[3 E3. 



fied Urea reactions is negligible compared to the photon 
emission from the surface of the star. However, if the 
star is massive enough to enable direct Urea reactions 
in the core (such as for our 2.21M neutron star model) 
then the neutrino emission will dominate the cooling pro- 
cess for surface temperatures higher than about 40 eV 
[38 . The thermal equilibrium condition can be written 
as Wd = L^ + Ly, where reheating due to r-mode dissipa- 
tion is given by Wd = an d is a function of r-modc 
amplitude, a. Therefore, a can be written in terms of 
luminosities as: 



5 x 3 4 



Q 



2 s JMR 3 fl i v 2ttG 



-) 



1/2 



(19) 



where L7 = 47ri? 2 crT 4 // is the thermal photon luminosity 
at the surface of the star. Here, R and T e ff are the 
stellar radius and surface temperature, respectively, and 
the neutrino luminosity is given by 

L = ^Rdu^qcdLdu t6 | MR 3 ~ Rdu)^QcdLmu t8 



A 4 



A 4 



where T is the core temperature, Rdu is the radius of 
the core where direct Urea neutrino emission is allowed 
and L is a dimensionless parameter given in Table |lj The 
thermal equilibrium condition for a NS with standard 
neutrino cooling (1.4M and 2.OM NSs in this study) 
can be approximated as Wd — since the neutrino 
cooling in this case is negligible compared to the surface 
photon luminosity. 

As can be seen in Figjl] , out of the 11 sources consid- 
ered in this paper only 9 of them are likely to have un- 
stable r-modes, meaning that they are above the r-mode 
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Source 


Ot-th.eq 

(1.4M a ) 


Ot-th.eq 

(2.OM ) 


Ot-th.eq 

(2.21M ) 


v (Hz s" 1 ) 
(1.4M ) 


v (Hz s^ 1 ) 
(2.OM ) 


v (Hz s^ 1 ) 
(2.21M ) 


v a d (Hz s 1 ) 
observation 


4U 1608-522 


7.15 x 10" 8 


6.60 x 10" 8 


1.18 x 10" 5 


-1.44 x 10" 15 


-1.78 x 10 -15 


-4.25 x 10" 11 




IGR J00291+5934 


1.14 x 10 -8 


1.32 x 10 -8 


1.82 x 10" 7 


-2.89 x 10" 17 


-5.59 x 10~ 17 


-7.94 x 10 -15 


-3 x 10 -15 


MXB 1659-29 


1.16 x 10~ 8 


1.07 x 10~ 8 


6.90 x 10~ 8 


-1.78 x 10~ 17 


-2.18 x 10~ 17 


-6.77 x 10~ 16 




Aql X-l 


3.49 x 10" 8 


3.27 x 10" 8 


1.02 x 10~ 6 


-1.49 x 10" 16 


-1.89 x 10 _le 


-1.37 x 10" 13 




KS 1731-260 


2.35 x 10" 8 


2.20 x 10" 8 


2.93 x 10~ 7 


-4.81 x 10" 17 


-6.09 x 10~ 17 


-8.07 x 10" 15 




XTE J1751-305 


5.09 x 10~ 8 


4.76 x 10~ 8 


6.54 x 10~ 7 


-6.13 x 10~ 17 


-7.74 x 10~ 17 


-1.09 x 10~ 14 


-5.5 x 10~ 15 


SAX J1808-3658 


1.28 x 10" 8 


1.19 x 10" 8 


2.16 x 10" 8 


-2.19 x 10" 18 


-2.74 x 10 -18 


-6.74 x 10~ 18 


-5.5 x 10" 16 


XTE J1814-338 


1.76 x 10~ 7 


1.67 x 10" 7 


2.04 x 10" 6 


-7.49 x 10" 17 


-9.73 x 10~ 17 


-1.08 x 10" 14 




NGC 6440 


1.54 x 10~ 6 


1.45 x 10 -6 


3.64 x 10~ 5 


-2.90 x 10~ 16 


-3.71 x 10~ 16 


-1.75 x 10~ 13 





Table IV: Upper bounds on the r-mode amplitude from the "thermal equilibrium" condition that are consistent with quiescent 
luminosity data are given for different star models. The gravitational radiation induced spin-down rates due to unstable r-modes 
as well as the observed spin-down rate for some of the sources are also given |16l 1301 155] . 



instability curve. The two most slowly rotating sources, 
XTE J1807 and XTE J0929, are outside the instability 
region for a hadronic star, which means they likely can no 
longer spin down due to gravitational radiation from an 
r-mode. Therefore, we only evaluate the upper bounds 
on the r-mode amplitude for those 9 sources within the 
instability window. Using Eq. |19| we have evaluated a 
for all of those sources using the three different NS mod- 
els considered in this work. The values of a are given in 
As can be seen for the 1.4M Q and 2.0M© NSs, 



Table IV 



where there is no enhanced neutrino emission, the values 
of a range from 1.14 x 10~ 8 to 1.54 x 1CT 6 , where NGC 
6440-which has the lowest spin frequency-has the highest 
r-mode amplitude, as expected. In the case of 2.21M 
NSs the upper bounds on a are larger since due to the 
direct Urea neutrino emission Wj, can be larger, and a 
in this case ranges from 2.16 x 10~ 8 to 3.64 x 10~ 5 . 4U 
1608-522 has the highest temperature among the sources 
considered which implies a very large neutrino luminosity 
for the 2.21M Q NS model and as a result a large r-mode 
amplitude. But it should be noted that this system has 
been accreting for a long time and both its high quies- 
cent luminosity and surface temperature may be due to 
the long term accretion and our assumption in ascribing 
all the quiescent luminosity to the heat that comes from 
inside of the star is likely not a good estimate for this sys- 
tem, but here we are only interested in obtaining upper 
limits on r-mode amplitude. 



C. r-mode spin-down 

To see whether these results are plausible or not and 
what fraction of the quiescent spin-down of these sources 
can be due to gravitational radiation from r-mode oscilla- 
tions, we use our results for a from the "thermal equilib- 
rium" condition and insert them into the right hand side 
of Eq. |12a| to determine v in the absence of accretion 



(N aC c — 0). The derived spin-down values for different 
neutron star models are given in Table [TV] and are shown 
graphically in Figure 2 (vertical green lines). Comparing 
these results with the observed spin-down rates, which 
exist for IGR J00291, XTE 1751-305 and SAX J1808 (red 
diamonds in Figure 2), we find that in the case of 1.4M 
("x" symbols) and 2.0M Q (square symbols) NSs, where 
there is no fast neutrino cooling present in the star, the 
r-mode spin-down can only provide about 1% of the ob- 
served spin-down rate, which means that other spin-down 
mechanisms such as magnetic-dipole radiation are re- 
sponsible for spinning down a 1.4M Q or 2.OM hadronic 
star with no fast cooling process. For the two AMXPs 
with relatively "slow" spin frequencies, NGC 6640 and 
XTE Jl 8 14-338, we find r-mode spin-down limits that 
are more competitive with observed values. Indeed, our 
limit for NGC 6640 is comparable to the measured, qui- 
escent spin-down rate for SAX J 1808, thus spin-down 
measurements for this source would be particularly in- 
teresting. 

On the other hand, the values of spin-down rate that 
we get for a 2.21Af (circle symbols) hadronic star with 
direct Urea neutrino emission is not consistent with the 
observations for most of the sources since it gives a very 
large rate and in the case of IGR J00291 and XTE J1814 
it is larger than the observed values. However, it is con- 
sistent for SAX J1808 and perhaps MXB 1659-29, and 
in the case of SAX J1808 r-mode spin-down is still about 
1% of the observed spin-down rate. 

With typical values of a few xlO -8 , our derived am- 
plitude upper limits suggest that for many LMXB neu- 
tron stars the r-modes are likely not excited to suffi- 
cient amplitudes to substantially affect their spin evo- 
lution. This begs the question of whether or not un- 
stable, steady-state r-modes actually exist in these neu- 
tron stars? One possibility is that additional damping 
mechanisms, such as perhaps associated with crust ef- 
fects, superfluid mutual friction [5] or the existence of 
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100 1000 

v (Hz) 

Figure 2: Limits on the spin-down rates due to an r-mode torque for nine LMXB systems and a range of neutron star masses 
are shown in the v versus v plane. The vertical green lines show the full range of v for each source (labelled) obtained from 
the r-mode amplitude limits derived from observations of quiescent luminosities and temperatures (see the discussion in §111) . 
The limits obtained for 1.4, 2.0, and 2.21M0 NS models are marked by the "a;," square, and circle symbols, respectively. Also 
shown are two pairs of parallel lines representing magnetic-dipole (solid) and r-mode (dashed) braking laws. Lines are drawn 
for two values of the magnetic field, 10 s G (lower), and 10 9 G (upper), as well as two values of a, 5 x 10~ 8 (lower) and 5 x 1CP 6 
(upper). For systems with measured, quiescent spin-down rates, these values are marked with the red diamond symbols. For 
additional context, millisecond pulsars from the ATNF pulsar database |40| are shown and denoted by the black "+" symbols. 



exotic matter in the core of NSs [23] , are at work and 
modify the instability window so as to render these neu- 
tron stars stable to r-mode excitation. An interesting, 
related question is whether the existence of r-modes at 
the amplitudes estimated here can be inferred directly 
from observations? Figure 2 shows both r-mode (dashed 
parallel lines) and magnetic-dipole (solid parallel lines) 
spin-down laws. The r-mode spin-down braking index, 
n = 7, is steeper compared to that for magnetic spin- 
down (n = 3) , thus at high enough spin frequencies (well 
above a kHz) one might expect that the r-mode torque 
would eventually become competitive with or dominate 
the magnetic dipole torque. However, as of yet there are 
no known neutron stars spinning fast enough for this ef- 
fect to become dominant, and depending on the EOS the 
mass-shedding limit might be reached before the r-mode 
torque becomes competitive with the magnetic torque. 
As discussed above, quiescent spin-down measure- 



ments have been typically attributed to the magnetic 
torque. For a number of sources considered here our r- 
mode amplitude limits support this presumption. Any 
spin-down contribution from an r-mode torque would be 
more easily identifiable if the magnetic field strengths of 
these neutron stars were constrained independently from 
the magnetic spin-down estimate. Moreover, identify- 
ing an r-mode spin-down would, in principle, be sim- 
pler for those neutron stars with the lowest magnetic 
torques, and thus field strengths. At present, the low- 
est inferred millisecond pulsar magnetic field strengths 
are w 6 x 10 7 G. At this level the magnetic spin-down 
is of order i) « 5 x 10 -17 Hz s _1 , which is compara- 
ble to our derived r-mode spin-down limits for a num- 
ber of sources considered here, assuming their neutron 
stars are < 2M . In this regard, quiescent spin-down 
measurements for more of the sources considered here, in 
particular the AMXPs XTE J1814 and NGC 6440, would 



9 



be extremely valuable. 



D. Gravitational wave amplitudes 

R-modes in NSs are one of the possible mechanisms for 
gravitational wave (GW) emission and they can be obser- 
vationally interesting in newborn neutron stars and per- 
haps accreting neutron stars in LMXBs [41 . Continuous 
GW emission from r-modes is dominated by I = m = 2 
current quadrupole emission [4 . The gravitational wave 
amplitude ho (strain tensor amplitude) is related to the 
r-mode amplitude a by the following equation |41l 142 j : 



worth mentioning that indirect upper limits on GW am- 
plitude can be obtained for sources with observed spin- 
down rates, i>, by assuming that all of the observed spin- 
down is due to GW emission |41j : 



h sd — 
n — 



1 45GIP 



8c 3 P 



(23) 



where P = is the observed pulse period and | p | = | ^ | . 
Using this equation for the three sources with measured v 
and different NS models (i.e. different masses and radii) 
we get hff" values that range from 4.14 x 10 -28 to 6.53 x 
lO" 28 . 



h = 



-aujlMR 3 ,}, 



(21) 



where r is the distance to the source, M and R are 
the mass and radius of the NS, J is the dimensionless 
parameter defined by Eq.|3]and ui r is the frequency of the 
r-mode, which is related to the angular spin frequency of 
the star f2 (for the m = I = 2 r-mode) by the following 
equation 



ui r 



-n 



(22) 



Using the upper limits on the r-mode amplitude of NSs 
in LMXBs derived above, we can obtain upper limits on 
the amplitude of the GWs emitted from these sources 
due to unstable r-modes. For the sources considered in 
this work, upper limits on the GW strain amplitude ho 
for the 1.4M and 2.OM NS models are in the range 
of 1.8 x 1CT 29 to 4.9 x 1(T 28 which is below the antici- 
pated detectability threshold of Advanced LIGO [3TJ [33] . 
In the case of the 2.21M© NS model, since the r-mode 
amplitudes are larger than those for the low mass NSs, 
we get larger values of ho, but even in this case for most 
of the sources ho is still below the detectability threshold 
of Advanced LIGO. The highest values of ho for a mas- 
sive star are obtained for NGC 6440 and Aql X-l with 
an amplitude of the order of 4 x 1CT 27 and 4U 1608-522 
with an amplitude of 7.18 x 10 -26 . However, it should be 
noted that the large values of r-mode amplitude in these 
sources, which would cause very large spin-down rates, 
may eventually be ruled out by future measurements of 
their spin-down rates. In this context it is important to 
restore the X-ray timing capability that was lost when 
RXTE was decommissioned in January 2012. Missions 
planned or currently in development which could pro- 
vide such a capability include India's ASTROS AT [44], 
ESA's Large Observatory for X-ray Timing (LOFT, |4"5]), 
and NASA's Neutron Star Interior Composition Explorer 
(NICER, 051). 

The upper limits on GW amplitudes discussed here 
are related to the GW emission due to unstable r-modes 
but our results do not exclude the possibility of having 
larger GW amplitudes in LMXBs from other GW emis- 
sion mechanisms such as NS mountains [4T1 US]. It is 



IV. CONCLUSIONS 

In this paper we have presented upper limits on r- 
mode amplitudes in LMXB neutron stars using their 
observed quiescent luminosities, temperatures and spin- 
down rates. We calculated results for NS models con- 
structed with the APR EOS (normal hadronic matter) 
with masses of 1.4, 2 and 2.21 M©, where our highest 
mass model (2.21 M©) can support enhanced, direct Urea 
neutrino emission in the core. We have used two differ- 
ent methods to calculate r-mode amplitudes. The first is 
based on the assumption that in an outburst-quiescence 
cycle all the spin-up torque due to accretion during the 
outburst is balanced by the r-mode spin-down torque 
due to gravitational radiation. This method gives am- 
plitudes in the range of « 1 - 3 X 10~ 7 for the sources 
with measured spin-up rates. Since in reality there are 
other sources of spin-down such as magnetic-dipole ra- 
diation that may be the dominant spin-down source, we 
use another method for computing the r-mode amplitude 
that doesn't ascribe all of the spin-down of the star to 
gravitational radiation and therefore gives tighter bounds 
on the amplitudes. This second method is based on the 
assumption that in a thermal steady state some fraction 
of the heat that is generated in the star due to r-mode 
dissipation will be lost from the star by neutrino emission 
and the rest will be radiated at the surface. This assumes 
that all of the heat emitted from the surface of the star 
during quiescence is due to the r-mode dissipation inside 
the star, and thus provides an upper bound on the r- 
mode amplitude. We have computed core temperatures 
as well as neutrino and thermal (photon) luminosities 
for LMXB sources using measurements of the quiescent 
luminosities and surface temperatures and showed that 
at temperatures relevant for LMXB neutron stars, when 
there is no enhanced cooling mechanism, the cooling of 
the star is dominated by photon emission from the sur- 
face (for 1.4 and 2.0M© NS models), but in a massive 
star where direct Urea neutrino emission is allowed, the 
cooling is dominated by neutrino emission (for T e f f > 40 
eV). For the lower mass NS models (1.4 and 2 M©) we 
find dimensionless r-mode amplitudes in the range from 
about 1 x 10~ 8 to 1.5 x 10~ 6 . For the AMXP sources 
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with known quiescent spin-down rates these limits sug- 
gest that < 1% of the observed rate can be due to an un- 
stable r-mode. Interestingly, the AMXP with the highest 
amplitude limit, NGC 6640, could have an r-mode spin- 
down rate comparable to the observed, quiescent rate for 
SAX J1808. Thus, quiescent spin-down measurements 
for this source would be particularly interesting. Having 
enhanced, direct Urea neutrino emission in the core of our 
highest mass model (2.21 M Q ) means that the dissipated 
heat in the star can be larger and therefore it can have 
higher r-mode amplitudes. Indeed, the inferred r-mode 
spin-down rates at these higher amplitudes are inconsis- 
tent with the observed spin-down rates for some of the 
LMXB sources, such as IGR J00291 and XTE J1751-305. 
This could be used to put an upper limit on the masses of 
these sources if they were made of normal nuclear matter, 
or it could be used to probe the existence of exotic mat- 
ter in them if the NS mass in these systems were known. 
Using the results for r-mode amplitudes, the upper lim- 
its on gravitational wave amplitude due to r-modes have 
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